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INTRODUCTION 


We  begin  with  a  theoretical  characterization  of  the  Cn  basis  functions 
with  which  we  are  dealing.  For  a  given  node  x-^  in  R,  there  is  a  basis 
function  H-y(x)  associated  with  the  jth  derivative  of  any  function  f  at  x^, 
where  j  ranges  from  o  to  n.  In  addition,  each  basis  function  is  nonzero  on 
only  two  adjacent  subintervals.  Continuing  with  the  definition  of  the  H's,  we 
wish  an  approximation  F  to  f  of  the  form: 

n  (j )  ( j ) 

F(x)  =  I  HijCxjf  (xi)  +  Hi+ij(x)f  (xi+i) 

j-o 

where  x^  <  x  <  x^-f  j . 

In  order  that  F  and  its  n  derivatives  will  agree  with  f  and  n  of  its 
derivatives  at  nodes  x^  and  x^+i ,  it  is  sufficient  that  the  basis  functions 
associated  with  arbitrary  node  i  obey  the  following  conditions: 

(k) 

Hij(xi)  =  Sjk 

and 

(k)  (k) 

Hij(xi-l)  =  Hij(xi+i)  =  0 

where  0  <  j  ,k  <  n. 

On  a  given  subinterval  each  H  must  obey  n+1  conditions  on  the  left 
extreme  and  n+1  conditions  on  the  right  extreme.  The  H's  can  therefore  be 
represented  by  two  distinct  polynomials  of  degree  2n+l . 

The  following  series  of  pictures  depicts  the  C1*  basis  functions  (scaled) 
and  their  derivatives.  The  five  functions  across  the  top  are  the  basis 
functions  associated  with  the  Oth  through  the  4th  derivatives  of  f,  and  the 
functions  underneath  them  are  their  successive  derivatives.  Note  that  the 


1 


functions  along  the  diagonal  are  nonzero  in  the  center  while  all  off  diagonal 


functions  are  zero  there. 


Figure  1.  C1*  Basis  Functions  and  Their  Derivatives. 


Although  subsequent  analysis  will  enable  us  to  compute  all  the  basis 
functions  for  any  given  level  of  smoothness  (n)  in  an  extremely  simple  way, 
the  author  has  not  seen  anything  similar  mentioned  or  referenced  in  any  finite 
element  text  thus  far. 

DERIVATION 

We  begin  by  defining  a  finite  support  Taylor  series  (FSTS).  Take  the 
ordinary  Taylor  series  for  f  around  node  x^,  truncate  it  beyond  nth  derivative 
terms,  and  multiply  each  term  by  a  function  which  will  have  the  effect  of  (1) 
not  disturbing  the  truncated  Taylor  series  at  all  at  node  x^  and  (2)  zeroing 
the  series  and  n  of  its  derivatives  at  nodes  x-^-q  and  x±+\ .  If  we  do  this  for 
each  node  ,  calling  the  result  F^(x) ,  we  get  a  global  approximation  to  f 
which  agrees  with  the  truncated  Taylor  series  of  f  at  each  node  by  simply 
summing  the  F^.  This  is  just  an  alternative  way  of  defining  the  piecewise 
Hermite  approximation  which  we  will  find  quite  useful. 

TS:  f ( x )  «  l  f  IxiKx-xjV /j  !  (D 

3  =  0 

n  (j )  j 

FSTS:  Fi(x)  »  l  f  (Xi)(x-Xi)  gj(Ri(x))/j!  (2) 

3=0 

where 

R^(x)  =  1  if  x  <  Xi-i  or  x  >  x^+i 

=  (xi-x)/(xi-xi-i)  if  xi-i  <  x  <  x± 

“  (x-xi)/(xi+1-xi)  if  Xi  <  x  <  Xi+i  (3) 

Rj^(x)  is  just  one  minus  the  hat  function  associated  with  node  i  or  just  the 
relative  position  of  x  in  either  the  left-  or  the  right-hand  subinterval.  The 
domain  of  the  g  functions  is  therefore  just  the  interval  [0,1]. 
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The  objective  now  is  to  determine  the  g's.  Since  we  want  F±  and  its 
derivatives  to  behave  in  a  certain  manner,  we  must  first  differentiate  F^x) 
an  arbitrary  number  of  times.  Using  Leibniz's  rule  for  differentiating  a 
product,  we  have: 

.  V  f(p  >  mlnh  J-k  <»-W  '  »-k 

Fx(x)  -  l  £  <X1)  l  (fcXx-x!)  gj(R1(x))(Ri(x))  /(j-k)!  (4) 

j=0  k=0 

and  substituting  x  =  x*  in  Eq.  (4)  we  have 

(m)  ®  ( j )  m  (m-j)  '  m-j 

Fi  (xi)  “If  (xi)(j)gj  (0X^(0)) 

j-o 

t  I 

We  can  define  R^O)  to  be  Ri(±e)  or  take  limits  from  either  side  of  Xi. 

We  now  want  conditions  on  the  g's  which  are  sufficient  for: 

(m)  (m) 

Fi(M)  =  f  (Xi)  (6) 

and 


(5) 


(m)  (m) 

Fi(xi-l)  =  Fi(xi+1)  =  0 


(7) 


for  0  <  m  <  n. 

We  can  get  these  conditions  from  Eqs.  (4)  and  (5).  Conditions  on  the 
g's  sufficient  for  Eqs.  (6)  and  (7)  are: 

gj(0)  =  1  0  <  j  <  n  (g) 


(m) 

gj  ( 0)  =  0  0  <  j  <  n,  1  <  m  <  n-j 

(m) 


(9) 


gj(l)  =  0  0<j<n,0<m<n  (10) 

The  2n-j+2  conditions  on  gj  can  be  met  by  a  polynomial  of  degree  2n-j+l. 
The  product  of  gj  and  (x-xi)j  in  Eq.  (2)  is  therefore  of  degree  2n+l  for  all 


j ,  as  expected. 


The  g  of  lowest  degree  is  gn,  which  has  defining  conditions: 


gn<0)  -  1 


and 


0  <  m  <  n 


This  g,  determined  by  inspection,  is: 


n+1 


gn(x)  “  (1-x) 


(11) 


Now,  from  Eq.  (10),  we  see  that  all  the  g's  have  the  same  derivative  behavior 
at  x=l.  We  then  need  only  define  gj (x)  as  the  product  of  some  unknown 
polynomial  hj (x)  and  gn(x) : 


n+1 

gj(x)  =  hj(x)d-x)  0  <  j  <  n 


(12) 


where  hj  is  a  polynomial  of  degree  n-j : 


n-j  ^ 

hj  (x)  =  l  a^x  0  <  j  <  n 

k=0 


(13) 


It  may  seem  at  first  glance  that  the  a’s  should  have  an  extra  subscript  - 
namely  j,  since  we  are  seeking  n+1  sets  of  coefficients.  As  will  become 
immediately  apparent,  however,  one  set  is  sufficient  and  all  other  sets  are 
subsets  of  this  one.  This  subset  property,  along  vrLth  the  fact  that  the 
largest  set  of  a's  can  be  obtained  in  an  almost  trivial  manner,  makes  the 
result  of  this  analysis  truly  simple  indeed. 

If  we  now  obtain  the  mth  derivative  of  gj (x)  and  evaluate  it  at  x  =  0,  we 


get: 


(14) 
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and 


Using  Eqs.  (8)  and  (9),  we  have: 

ao  38  1 

l  ak(-l)m”k(n+1)  =  0  for  1  <  m  <  n-j  (15) 

k=0  m-k 

This  is  a  lower  triangular  system,  which  can  be  easily  solved  for  the  a's  by 
forward  substitution.  Note  that  the  coefficients  do  not  depend  on  j;  so  we 
might  as  well  solve  the  largest  system  (j=0)  and  obtain  all  the  a's,  although 

only  the  first  n-j+1  a's  are  needed  for  hj(x). 

Solving  system  (15)  for  j=0  and  a  few  values  of  n  gives  us  the  following 

table  of  afs 

TABLE  I.  COEFFICIENTS  OF 


Inspection  of  this  table  gives  us  a  very  simple  recursion  for  the  n  set  of 
coefficients  in  terms  of  the  n-1  set: 

n 

a0  =  1 

n  n  n-1 

ak  13  ak— 1  +  ak  1  ^  k  <  n— 1 

n  n  . 

an  =  2an_i  (16) 

where  the  superscripts  denote  the  level  of  smoothness* 


6 


We  therefore  have  here  nothing  more  than  one-half  of  Pascal's  triangle, 
viewed  at  an  angle! 

Recalling  that: 

n+1 

gj(x)  =  hj(x)(l-x) 

and  from  the  FSTS  that: 

Hij(x)  =  gj  (Ri(x)  )(x-x^)^ /j  ! 


we  have,  explicitly: 

n+1  j 

Hij(x)  =  hj  (Ri(x)  )(  l-R^(x)  )  (x-x^)  /j  ! 

Therefore,  f  may  be  approximated  on  [x-j^xi+i]  by: 

n  ( j )  x-xi  xn-j-x  n+1 


5  f  (j)  x-xi  xi+1-x  n+1  j 

1  {f  (xi)hj  ( - )( - )  (x-x^) 

j=o  . - . . 


xi+i-xi  xi+1-xi 
xi+l~x  x-xi  n+1 


Ci;  xi+i-x  x-xi  n+i  j 

+  f  Cxi+1)hj( - )( - )  (x-xi+1)  }/j  ! 

Xi+l-Xi  xi+1-xi 


(17) 


(18) 


CONCLaSlON 

In  order  to  evaluate  the  derivatives  of  the  H's,  one  can  expand  the 
polynomials  involved  and  multiply  out  or  one  can  apply  Leibniz’s  rule  a  couple 
of  times.  The  latter  course  is  deemed  simpler  and  more  numerically  stable 
since  it  leaves  a  result  which  is  in  "nearly"  fully  factored  form.  The  latter 
method  was  used  to  produce  Figure  1. 
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